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Abstract. This paper is concerned with the stabihty of steady sohitions to initial- 
boundary-value problems of reaction-hyperbolic systems for axonal transport. Under 
proper structural assumptions, we clarify the relaxation structure of the reaction-hyperbolic 
systems and show the time-asymptotic stability of steady solutions or relaxation boundary- 
layers. 
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1. Introduction 
This work is concerned with the following reaction-hyperbolic system 

r 

(1.1) dtUi + Xida:Ui = ^kijUj, z = l,2,---,r 

i=i 

on the quarter-plane x,t > 0. Here Ui = Ui{x,t){i = 1,2, ■■ ■ , r) are unknown functions, 
Aj(z = 1, 2, ■ ■ ■ , r) and kij{i,j = 1, 2, ■ ■ ■ , r) are given constants. It is well believed [TU] 
that such models describe the axonal transport in neuroscience. 

The axonal transport is important for the maintenance and functions of nerve cells or 
neurons. A neuron consists of three parts mainly: cell body, dendrites and a single axon. 
The axon is a long and thin pipe whose length can exceed 10,000 times its diameter. It 
is this axon that distinguishes neurons from cells in other organs or tissues. The axon is 
responsible for signal transmission in the nervous system. Its cytoplasm does not contain 
rough endoplasmic reticulum and therefore its proteins can only be transported from the 
cell body, where all proteins are synthesized. In addition, the axonal transport is also 
used by the neuroscientists to trace the connections in the brain. 

The transport proceeds as follows. Proteins are stored in vesicles as cargos. The vesi- 
cles are attached to kinesin (anterograde motors) or dynein (retrograde motors) proteins. 
These motor proteins drive the vesicles to walk along the cytoskeletal microtubules as 
track. Here the kinesin proteins move the vesicles from the cell body to synapse (antero- 
grade transport), while the dynein proteins move the vesicles in the opposite direction 
(retrograde transport). During the transport, many biochemical processes are possible. 
For example, the cargos can leave its track, can switch its motor proteins from kinesin 
to dynein or vice verse, and can move back onto the track. Thus, we can divide the 
cargos into a number of subpopulations, such as free vesicles, vesicle-kinesin compounds 
off track, moving vesicle-dynein compounds on track, etc. 

As the axon is long and thin, it is reasonable to assume the transport only along the 
longitudinal direction of the axon. Denote by x > the distance down the axon from the 
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cell body which is located at x = 0. Let Ui = Ui{x, t) be the concentration at space-time 
{x,t) of the i-th subpopulations. According to Reed and Blum |10], the mathematical 
model for axonal transport is partial differential equations of the form 

(1.2) dtUi + Xid^Ui = Fi{ui,U2,- ■ ■ ,Ur), z = l,2,---,r 

defined on the quarter-plane x,t > 0. Here the term XidxUi accounts for the transport 
of the 2-th subpopulation with constant velocity Aj, and Fi{ui,U2, ■ ■ ■ ,Ur) describes the 
biochemical processes of the constituents. It is well recognized that the biochemical 
processes are much faster than the transport in biosystems. Thus, it is more proper to 
rewrite the general equation (11.21) as 

(1.3) dtUi + Xid^Ui = ^Fi{ui,U2,- ■ ■ ,Ur), z = l,2,---,r. 

Here the small parameter e > characterizes the fact that the biochemical processes are 
much faster than the transport. 

On the other hand, in a typical experiment for investigating axonal transport, a large 
quantity of radiolabeled amino acids are injected into a nerve ganglion. Thus, the proteins 
are synthesized with the injected amino acids in the cell body and transported along the 
axon, and finally the radioactivity appears in the axon in a few hours. The wave of 
radiolabeled proteins travels at approximately constant velocity. The shapes and speeds 
of the profiles depend on nerves, animals, and proteins being studied. Usually the amino 
acids are continuously available in the experiment, and the sharp approximately parallel 
the wave fronts which suggest traveling waves. There are at least two such systems which 
carry membrane-bounded organelles and cytoplasm separately. For more details about 
the experiments for axonal transport, the reader can refer to [H [3, [5]. 

In order to explain the approximate traveling waves observed in experiments, the au- 
thors of [21 El El E] studied the linear case (II. ip . where Fi{ui,U2, ■ ■ ■ ,Ur) is linear with 
respect to the Uj^s. A concrete example can be found in [3J on neurofilament transport. 
Because parabolic equations typically admit traveling wave solutions, these authors all 
related the system (11.31) to a parabolic equation — the diffusive limit of (II. 3p . Especially, 
in [6] Friedman and Hu used parabolic-type estimates to analyse the diffusive limit of the 
linear systems. However, it is known [Hill] that, unlike the conservation laws, hyperbolic 
systems with relaxation also allow non-trivial traveling wave solutions. 

In this project, we intend to explain the experimental observation by directly studying 
steady solutions to the initial-boundary-value problems of the above hyperbolic-reaction 
systems. Steady solutions very much look like the traveling waves but they should called 
(relaxation) boundary-layers due to the presence of the boundary x = 0. Because they can 
be observed experimentally, we believe that the steady solutions are stable. Therefore, the 
goal of this paper is to investigate the time-asymptotic stability of the steady solutions 
under proper structural assumptions on the system. 

The standard structural assumptions on the axonal transport models (II. ip read as 

(HI). %>0if2^j; 

(H2). (conservation of mass) J2l=i ^ij = i = 2, ■ • • , r; 

(H3). (irreducibility) for any i\, there is a sequence of indices ji, j2, ■ ■ ■ -.31 
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such that io = ji, ii = ji and kj^j^^^ > for 1 < m < / — 1; 

(H4). there exist i and j such that Aj ^ Xj. 

These assumptions are taken from [6]. See also [21 El El [lOl [H]- Remark that (H4) 
ensures the hyperbohc-reaction system fll.ip won't degenerate into a system of ordinary 
differential equations. 

As a first step of this project, we consider the linear system (II .ip and additionally make 
the following technical assumption in this paper that 

(H5). A, >0, z = l,2,---,r. 

In the future, we will remove or relax this (H5). 

With the above assumptions, we formulate the problem as follows. Write 

U = {ui,U2,- ■ ■ ,Ur)^, A = diag(Ai, A2, ■ ■ ■ , Ar), K = [kij]nxn- 

Then fll.ip can be written as 

(1.4) Ut + AU, = KU. 
The initial and boundary data are 

(1.5) U\t=o = Uo{x), t/U=o = f/o(0). 

Here the second equality implies the consistency condition which is necessary for the above 
initial-boundary- value problems to have a continuous solution. Moreover, we assume that 

(1.6) AUoM = KUo{0), 

which is necessary for the above problems to have a continuously differentiable solution. 
Our main results are the time-asymptotic stability of steady solutions to (11. 4p together 
with f lL5|) . 

For references on nonlinear systems for axonal transport, we mention that Carr in 
[2] discussed the existence of global classical solutions to a class of nonlinear models. 
Recently, in [T7] we clarified the relaxation structure of nonlinear models in [2] and verified 
the relaxation limit of BV-solutions to the Cauchy problems. 

The rest of this paper is organized as follows. In Section 2 we discuss the relaxation 
structure of the hyperbolic-reaction system (II. ip . Section 3 is devoted to the existence 
and stability of the steady solutions for (II. 4p . 

2. Relaxation Structure 

In this section, we show that the hyperbolic-reaction system (II. ip possesses the relax- 
ation structure formulated in provided that the assumptions (H1)-(H3) hold. We 
start with the following elementary fact as (a) of Lemma 1 proved in [11]. 

Lemma 2.1. Under the assumptions (H1)-(H3), the kernel of the matrix K is one- 
dimensional and is spanned by a vector with strictly positive entries. 

On the basis of this fact, we can show 

Lemma 2.2. is a single eigenvalue of K. 
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Proof. Set 



1 ei 

Ir-l 



where Ci = (1, 1, ■ ■ ■ , 1) is a vector in W~^, I^-i is the unit matrix of dimension (r — 1), 
and the superscript T denotes the transpose of vectors or matrices. Then the inverse of 
Li is 



-1 



1 -ei 

Ir-l 



Using the partition of Li, we rewrite K as 



K 



a 

/3 



where a = (fci2, fcis, ■ ■ ■ , kir), {3 = {k2i, k^i, - ■ ■ , Ki)'^ , and Ki is the (r— 1) x (r — l)-matrix. 
By a direct calculation using the assumption (H2), we obtain 



LiKLl 



-1 





/3 K2 



1 ei \ f kii a \ f I -ei 
^ Ir-l )\ 13 ^1 ; V ^r-i 
with K2 = Ki — [iei. 

Thus, it reduces to show that the sub-matrix K2 is invertible. Otherwise, there is a 
vector Tj = {r)i,r]2, - ■ ■ , ?7r-i)^ such that K2'r] = 0. Take 





-1 



V 



Then it holds that 



LI 



-1 



LI 



-1 







0. 



W 
P K2 J[v 

This shows that ^ is a vector in the kernel of K. On the other hand, we have 

/ \ ( -r]i-r]2 rjr-i \ 



i = L-i' 



1 -ei 

Ir-l 



V2 



\ 



Vr-l 



J 



\ Vr-1 f 

Obviously, such a ^ can not be in the kernel of K spanned by a vector with positive entries. 
This contradicts Lemma 12.11 Therefore, K2 is invertible and is a single eigenvalue of 
K. This completes the proof. □ 



Remark 2.1. Following the above proof, we have 







K2^fi /, 



r-l 



LiK 




K2 



1 

K2^ii Ir-l 



Li. 



On the other hand, it is not difficult to deduce from the Gershgorin circle theorem that 
non-zero eigenvalues of K have negative real parts. Therefore, K2 is stable and (i) of the 
first stability condition in [TB] is verified. 



Furthermore, we have 
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Lemma 2.3. Under the assumptions (H1)-(H3), there exist a positive definite diagonal 
matrix Aq, an orthogonal matrix P, and a positive definite diagonal matrix S such that 



AoA = AAo 
A,K + K^A, = -P^ 





s 



Furthermore, the first column of is in the kernel of K. 



P. 



Proof. Let ^ = {CiiC2,--- ,^r)^ is an eigenvector of the matrix K, associated with the 
eigenvalue 0. By Lemma [2. 1^ we may assume that > for each i. Define 



D 



6 



V J 

It is obvious that the matrix KD satisfies the assumptions (H1)-(H3) as well and 

= K 



KD 



6 



0. 



Namely, the sum of each row of KD is also 0. Moreover 
symmetric matrix KD + DK^ fulfils the (H1)-(H3), too. 

Take = D^^ and it is clear that Aq is positive definite diagonal and AqK 
Moreover, the symmetric matrix 



it is easy to see that the 



AqK + K^Aq = D-\KD 



DK^)D- 



also fulfils the assumptions (Hl)-H(3). According to Lemmas 12.11 and \2.2\ is a single 
eigenvalue of AqK+K'^Aq. Moreover, from the Gershgorin circle theorem it is not difficult 
to deduce that non-zero eigenvalues of AqK + K'^Aq are negative. Thus, there exist an 
orthogonal matrix P and a positive definite diagonal matrix S such that 



AoK + K^Ao = -P^ 





S 



P. 



Furthermore, since the r-vector (^1,^2, ■ ■ ■ ,^r)'^ is in the kernel of AqK + K^Aq, it is 
easy to see that the last (r — 1) components of the column vector F(^i,^2, ■ ■ ■ ,^r)"^ are 
zeros. Namely, the vector (^1,^2,- ■■ ,Cr) is orthogonal to the last (r — 1) rows of the 
orthogonal matrix P and thereby parallels to the first row of P. Hence, the first column 
of P^ is in the kernel of K and the proof is complete. □ 



Remark 2.2. Lemma [2.3[ together with Remark 12. H shows that the reaction-hyperbolic 
systems satisfying assumptions (H1)-(H3) fulfil the first stability condition in [T5] . 
However, they do not satisfy the second stability condition in [TJl |Tn] in general, unless 
further assumptions are posed. An important case is that KD is symmetric, which implies 
the second stability condition due to Theorem 5.3 in [15]. It is clear that KD is symmetric, 
provided that the principle of detailed balance holds (see, e.g., [13]). For the neurofilament 
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model in [3j, K is tri-diagonal and one can easily see that KD is symmetric. However, 
the assumptions (H1)-(H3) do not imply the symmetry of KD. In fact, the 4 x 4-matrix 

/ -4 1 1 \ 
2-301 
2 1-20 
V 1 1 -1 / 

satisfies (H1)-(H3) and (|, 1, 1, 2)^ is an eigenvector associated with 0. Set D = diag{^, 1, 1, 2). 
By a direct calculation, we have 



K 



KD 



( -2 
1 
1 

\ 



1 




\ 

2 


-2/ 



which is not symmetric. 
But we have 

Proposition 2.4. For r < 3, the matrix KD is symmetric. 
Proof. For r = 2, let X be 

K = 

with a,b > 0. We take ^ = (6, a)^ and D = diagip, a). Then 

KD = 



ab —ah 
—ah ab 



is symmetric. 
For r = 3, set 



K 



11 


kl2 


ki3 


21 




k23 


31 


h2 


kss 



and D = diag(^i, ^3)- Moreover, set 

a = 6^12 - 6^21, b = isku 

Since 



6^31, c = 6A;23 - 



32- 




6^11 + 6^12 + 6^13 
6^21 + 6^22 + 6^23 
6^31 + 6^32+6^33 

we see immediately from the assumption (H2) that 

a + b = b + c = c + a = 0, 

and therefore 

a = b = c = 0. 



0, 
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Hence the matrix 

KD = 

is symmetric and the proof is completed. □ 

Finally, we conclude this section with the following Lemma. 

Lemma 2.5. Under the assumptions (H1)-(H4), there is a skew symmetric matrix H 
and a positive constant c such that 

Proof, let ^ = (^1, ^2, • ■ ■ ) ?r)^ be in the kernel of the matrix K and Ao be the positive 
definite diagonal matrix in Lemma 12.31 Consider the system of equations 

AoWt + AqAW^ - {AoK + K^Ao)W = 0. 

From the proof of Lemma 12. 3[ we know that the kernel of the symmetric matrix AqK + 
K^Aq is equal to that of K. By Lemma 12. we may assume that > for each i. 
Thanks to the assumption (H4), one cannot find any number A such that 

for all i. Namely, ^ = (1^1,1^2, ■ ■ ■ ,^r)'^ is not an eigenvector of the coefficient matrix A. 
According to Shizuta and Kawashima (see Theorem 1.1 in [12] and also Theorem 2.3 in 
[16]). there exists a positive constant c and a skew symmetric matrix H such that 

HA- AH > cl + AoK + K^Aq. 

Hence the lemma is proved by combining this with Lemma 12.31 □ 

3. Existence and Stability of steady solutions 

In this section, we discuss the existence and stabihty of steady solutions of (11.41) . The 
equations for steady solutions B = B{x) are 

(3.1) AB^ = KB. 
Since the matrix A is invertible, W = AB satisfies 

(3.2) = {KA-^)W. 

From the assumption (H5) and the definition of A, we see that the matrix KA~^ satisfies 
the assumptions (H1)-(H3). Thus, we can use Lemma 12.21 and the Gershgorin circle 
theorem to show that KA~^ has (r — 1) stable eigenvalues and a zero-eigenvalue. 
Given any boundary data f/o(0), the solution for 03.21) is 

W{x) = exp{KA-^x)AUo{0), 

and therefore, 

B{x) = A-^ exp{KA-^x)AUo{0). 
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Since KA.~^ has (r — 1) stable eigenvalues and a zero-eigenvalue, the matrix exp(irA~^x) 
is bounded with respect to x > 0. In conclusion, given any boundary data -B(O), the 
equations (13. ip for steady solutions have a unique bounded solution B = B{x). 

Now we turn to discuss the stability. Set = U{x,t) — B{x). From (11. 4p and 

(13. ip . we see that $ satisfies 

(3.3) $t + A$^ = K$. 
The initial and boundary data are 

^{x,0)= Uo{x)-B{x), 

(3.4) $(0,t)= 0. 

From f lL6|) and i^J^, it follows that 

A{UoM - BM) = K{Uo{0) - B{0)). 

Since Uo{0) = B{0) and A is invertible, the consistency of the initial and boundary data 
becomes 

(3.5) $,,(0,0) = $(0,0) = 0. 

Thus, our task is reduced to analyzing time-asymptotic behaviors of the solution $ to the 
IB VP (Q together with 1^. 

We start with the following local existence result. This result can be showed by slightly 
modifying the proof given in Section 5 of [8] and we omit it here. 

Lemma 3.1. Suppose Uq{x) — B{x) G and i\3. 5\) holds. Then there exist a positive 
constant such that liS. 3\) together with Iji3.4\ ) has a unique solution $(x, t) G C(0, T^,; H"^). 
Moreover, the solution satisfies the following estimate 

sup ||<l>(-,t)||^2<2||[/o(x)-S(x)||^,2. 

0<t<T, 

Here depends only on the range of B{x) and any upper bound of \ \Uo{x) — B{x)\\h^ . 

Theorem 3.2. Under the assumptions (H1)-(H5), if Uo{x) — B{x) G and U.6\) holds, 
then together with U.5\) has a unique global solution U G C{0,oc; H"^) satisfying 

lim sup \U{x,t) - B{x)\ = 0. 

Proof. Let be the matrix given in Lemma 12.31 We multiply (13. 3 p with ^^Ao to get 

^'^Ao^t + ^'^AqA^,^ = ^'^AqK^. 

Namely, 

$^Ao$t + $'^v4oA$^ = ^<i>'^{AoK + K^Ao)<^> < -^'^^P' 
where we have used Lemma [2731 Setting 

V = P^ = {Vi,V2,--- ,Vrf 
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and integrating the above inequality with respect to (x, t) G [0, oo) x [0, t], we get 

+00 r+oo ^ 

-^^Ao^{x,t)dx - -^^AQ^{x,0)dx 
2 Jq 2 

/"* 1 c /"* 1'+°° ^ 

- AoA$(0, r)dr+ - / / y2Vi^{x,T)dxdT < 0. 

Jo ^ 2 Jq Jq .^^ 

Since Aq is positive definite and $(0,t) = 0, there is a generic constant C such that 

(3.6) + E fm{;T)\\,.'dr<Cm;0)\\l^- 

i=2 

Next we estimate the derivatives. Differentiating the equation (13. 3p with respect to x 
gives 

(3.7) ^xxt ~l~ -^^xxx -^^xx 

From the equations (13. Sp and (13. 7p with the data in ( 13. 4p . the boundary data for the 
derivative are 

$,(0, t) = A-\K^O, t) - $^(0, t)) = 0, 
$,,(0,t) = A-i(K<l>,(0,t)-$,i(0,t)) =0. 

Thus, we use the same technique shown above to estimate ||K'('i^)|| and ||14a;(',i^)||5 
ll<^>.(-,^)lli. + E f\\VU;r)\\L^"dT<C\\^,i.M\h, 

i=2 

W^U-Ml^ + Y. / ll^-x(-'^)IU^'^^<^ll*-(-'0)||i2. 

Summing these and the inequahty in (13.61) . we get 

(3.8) m-Mm + E r \m:r)\\lAT < C||<|.(-,0)||^.. 

i=2 

On the other hand, we multiply the equation (13. 3p with to get 

(3.9) ^iH^t + C^A$^ = ^iHK^. 
Here H is the skew symmetric matrix in Lemma 12.51 Since 

we integrate (13. 9p with respect to (x,t) and use Lemma [2751 to obtain 

ft p + OO ft f + OO 



c / / ^l{x,T)dxdT < / / Vji{x, T)dxdT 

Jo Jo Jo Jo 

+2 ^lHK^x,T)dxdT - / <^lH^{x,t)dx+ / <^lH^{x,0)dx. 

Jo Jo Jo Jo 
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Recall that $ = P^V thanks to the orthogonality of P and the first column of P^ is in the 
kernel of K. Thus, = KP^V is independent of the first component of V . Therefore 
the last inequality becomes 

ft p+oo pt p + oo ''' 

c / ^l{x,T)dxdT < / / 'S^V-^{x,T)dxdT 

Jo Jo Jo Jo 

pt p+oo 

+2 / / ^lHKP^V{x,T)dxdT 
Jo Jo 

p+oo 

^lH^{x,t)dx+ / (^lH^{x,0)dx 



< 



^^H^{x,t)dx+ I ^li 
Jo 

^ f\\^^{.,T)fdr + C r^r.(-,r)|||.dr 

^ ^0 Jo 



Namely, 

ft p+oo pt r 



pt p+oo Pt ' 

/ / ^l{x,r)dxdT<C $^||^.(-,r)||^.dr + C||$(-,t)||^. + C||$(-,0)||^.. 

Jo Jo Jo 

Combining this with (]3.8p gives 

(3.10) m-Mm + fY.\m:r)\\l,dr+ f\m-,r)\\Ur<Cm.M\l^ 

Jo "^0 

for all t > 0. 

Furthermore, for any ti,t2 > 0, we have 

\\\^x{-,t2)\\L^ - < \\M-,t2) - M-,tl)\\L^ 

rt2 pt2 

= 11/ A$^.^.(-,r)dr- / K^^{-,T)dT\\L2 
Jti Jti 

< C|t2 - til max ||$(-,t) 11^2 < C\t2-h\, 

T 

where we have used (13.71) in the second step. Thus, we can deduce that ||$^(-, t) ||£^2 — )■ 
as t — )■ +00. Moreover, from (13.101) . we see that ||$||l2 is bounded. Following from the 
celebrated Gagliardo-Nirenberg inequality 

|$|oo < V^||$||i||$.||i, 

we derive |$(-,t)|oo — )■ as t — t- +oo. Recall that $(x,t) = U{x,t) — B{x), and we 
complete the proof. □ 
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